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I. INTRODUCTION 

The description of nuclear structure and reactions in 
the unified theoretical framework is the long-standing 
challenge of nuclear theory. The attempts to reconcile the 
shell model (SM) with the reaction theory [l|, inspired 
the development of the continuum shell model (CSM) Q 
which evolved into the unified of theory of nuclear struc¬ 
ture and reactions [m. 

Structure of weakly bound states and resonances is 
different from the well-bound states. A comprehensive 
description of these systems goes beyond standard con¬ 
figuration interaction model such as the SM and requires 
an open quantum system formulation of the many-body 
system. Such a generalization of the standard SM to 
describe well bound, weakly bound and unbound many- 
body states is provided by the Gamow shell model (GSM) 
[§-[lOl|. GSM offers the most general treatment of cou¬ 
plings between discrete and scattering states. The many- 
body states in GSM are given by the linear combination 
of Slater determinants defined in the Berggren ensemble 
of single particle states which consists of Gamow (reso¬ 
nant) states and the non-resonant continuum. 

In this formulation, GSM is the tool par excellence for 
studies of the structure of bound and unbound many- 
body states and their decays. For the description of re¬ 
actions, the GSM has to be formulated in the CC repre¬ 
sentation. Recently, the GSM-GG has been applied for 
the calculation of excited states of ^®Ne and ^®Na, excita¬ 
tion function and the elastic/inelastic differential cross- 
sections in the ^®Ne(p,p') reaction at different energies 
[ni[ii. In this work, we apply the GSM-CG formal¬ 
ism for the description of low-energy radiative capture 
reactions: ^Be(p, 7 )®B and ^Li(p, 7 )®Li. In light nuclei, 
GSM-CG can be applied also for the description of nu¬ 
clear reactions in the ab initio framework of the no-core 


GSM [T^ and to heavier projectiles like the a-particle. 

The solution of solar neutrino problem is passing 
through an understanding of the ^Be(p, 7 )®B proton cap¬ 
ture reaction. ®B produced in the solar interior in this 
reaction, is the principal source of high energy neutri¬ 
nos detected in solar neutrino experiments. At the solar 
energies ( ~ 20 keV), this cross-section is too small to 
be directly measurable. For this reason, the theoretical 
analysis of this reaction is so important. On the other 
hand, whenever the measurement is feasible (Acm >150 
keV), the exact value of the capture cross section de¬ 
pends: (i) on the normalization obtained indirectly from 
the ’’Li(d,p)®Li cross section and, (ii) on the model de¬ 
pendent extrapolation of measured values of the cross- 
section down to the interesting domain of solar energies. 

Proton radiative capture reaction on ^Be is of particu¬ 
lar importance in astrophysics since it is involved in the 
pp-II and pp-III reaction chains. Indeed, the relative 
rates of the ^Be(e”, z/e)^Li reaction and the ^Be(p, 7 )®B 
reaction determine the pp-I/pp-II branching ratio, and 
thus the ratio of the neutrino fluxes coming from ^Be 
and ®B [I^. ^Be(p, 7 )®B reaction has been studied ex¬ 
perimentally by the direct proton c^ture [Tsl - I^ and 
the Coulomb dissociation of ®B [261 - K^ . Theoretical ap¬ 
proaches used to describe this reaction include the poten¬ 
tial model [Hj , the R-matrix approach [H, [33. the shell 
model embedded in the continuum (SMEC) [34| . the mi¬ 
croscopic cluster model [3^ , and the approach combining 
the resonating-group method and the no-core shell model 

M- 

^Li(n, 7 )®Li reaction is the mirror reaction of 
^Be(p, 7 )®B. ^Li(n, 7 )®Li reaction cross section 

at very low energies provides the essential ele¬ 
ment of rapid process of primordial nucleosynthe¬ 
sis of nuclei with A > 12 in the inhomogeneous 
big-bang models (37l-l40l|. Indeed, in the inhomoge- 
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neous big-bang hypothesis, the main reaction chain 
leading to the synthesis of heavy elements is [s^ 
^H(n, 7)-^^H(n, n)->^He(t, 7)-*-^Li(n, 7)®Li, 

and then ®Li(a, 

..., etc., for heavier nuclei. In this sense, the reaction 
^Li(n, 7 )®Li is a key process to bridge the gap of mass 
A = 8 and to produce heavy elements. The reaction 
^Li(n, 7 )®Li has been studied experimentally (dll - l^ . 
Theoretical studies of this reaction has been done using 
various potential models 1^ . the SMEC [s^, the 
microscopic cluster model and the halo effective 
field theory approach (d^ . 

The paper is organized as follows. The chapter [IT] 
presents the general formalism of the GSM-CC approach. 
In Sec. Ill Al we introduce the translationally-invariant 
GSM Hamiltonian in the cluster-orbital shell model 
(GOSM) variables [d^. The coupled-channel equations 
of the GSM-GG are presented in Sec. Ill Bl The channel 
states expansion in Berggren basis and the calculation of 
Hamiltonian matrix elements are discussed in Sects. EC] 
and III PI respectively. In Sec. Ill El we discuss how to 
orthogonalize the channel states, and Sec. Ill FI presents 
the method of solving the GSM-CC equations derived in 
this work. 

The chapter ED is devoted to the presentation of nu¬ 
cleon radiative capture formalism in the context of GSM- 
CC. In particular, the method of calculating matrix el¬ 
ements of the electromagnetic operators is explained in 
Sec. IIII Al and the matrix elements itself are given in the 
appendix. 

Results of the GSM-CC calculations are discussed in 
the chapters IIVI and IVl for ’’Be(p, 7 )®B and ’’Li(p, 7 )®Li 
low-energy reactions, respectively. Finally, main conclu¬ 
sions of this work are summarized in chapter I VII 


II. COUPLED-CHANNEL FORMULATION OF 
THE GAMOW SHELL MODEL 


A. Hamiltonian of the Gamow Shell Model 


spin-orbit term The Coulomb field is generated by a 
Gaussian density of Zc protons of the core i- vih-fj) 
in ([T|) is the two-body interaction which splits into nu¬ 
clear and Coulomb parts. As in the standard SM, adding 
and substracting a one-body mean-field U{fi) to the core 
Hamiltonian and the two-body interaction, respectively, 
allows to recast the GSM Hamiltonian in the form: 

H = Hbasis + T + Vj-es (2) 

where the potential C/basis generates the s.p. basis, the 
kinetic term is written T and the residual interaction is 
given by V^es- 


Ub 


asis 


A''val 


E (Ucifi) +U(fi)) 


( 3 ) 



In the present studies, we use the Furutani-Horiuchi- 
Tamagaki (FHT) finite-range two-body interaction [H, 

M- 


V{f 




( 5 ) 


The central potential is: 


V- = y In 

Z-/ 0, 






-H^P^ 


-M^P-P^) , 


( 6 ) 


where r is the distance between particles i and j, /3® is 
the range of gaussians, P'^ and P'^ are the spin exchange 
and isospin exchange operators, respectively, and 
By and are the exchange parameters. The spin- 
orbit potential writes: 


Center-of mass (CM) excitations in SM wave func¬ 
tions are removed using Lawson method (sol - l^ . In the 
Gamow shell model (GSM), this method cannot be used 
because Berggren states are not eigenstates of the har¬ 
monic oscillator (HO) potential. To avoid spurious CM 
excitations in GSM wave functions, the GSM Hamilto¬ 
nian is expressed in the intrinsic nucleon-core coordinates 
of the GOSM [H: 


Afval / ^2 


H = 


\ Wval / 

.)bE 

/ ^<J \ 


yin-y) + 


P^-Pj 

M, 


( 1 ) 


where Aval is the number of valence nucleons, Me is the 
core mass, and: llp.i = llMc+ ^/rrii, is the reduced mass 
of the i-th nucleon. The single-particle potential C4(f) 
which describes the field of the core acting on each nu¬ 
cleon, is a sum of nuclear and Coulomb terms. The nu¬ 
clear term is given by a Woods-Saxon (WS) field with a 


V.SO 


L-S y 

n=l 


SOg-/: 


■(IF; 


so 


pSO pr 


( 7 ) 


where L is the relative orbital angular momentum be¬ 
tween the two particles and S = Si + sj where Si , Sj are 
the spins of particles i, j. The tensor potential IdJ writes: 


Vj = Ot E (IF„^ - h:^P1 

n=l 


( 8 ) 


where 


Ot = 




- CTi-cr, r 


(9) 


and (Ti, aj are the Pauli matrices. The Coulomb potential 
in the FHT interaction is standard. 
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It is convenient to rewrite the FHT interaction using 
projection operators on singlet and triplet states of spin 
and isospin: 


|r, c) = yl(|r) ® |c)). In the above equation, (rjcldt) 
are the antisymmetrized channel wave functions: 
'l'c(r) E (r, cldt) E Uc{r)r. Hence: 


V^(t) = vS/C(r)« + vC/C(r)« 


( 10 ) 


ySO(r) = (11) 


,SO fSO 


so pSO, 


1^) = ^^ |r,c) . (15) 

C 

GSM-CC equations are obtained by inserting (HSJ in 
the Schrddinger equation and then projecting this equa¬ 
tion on a given channel basis state {r',c'\. One obtains: 


V^\r) = Ot(V,TaT(^)« + V-At (r)«) (12) 


T pT/ 


^ drr^ {Hc',c{r',r) - ENc',c{r',r)) =0 (16) 


where are the projection operators on singlet and 

triplet states of spin: 


where: 

Hc'Ar',r) = {r',c\H\r,c) 


(17) 


TTg ^ Si'Sj , 


Tt = -+s,-Sj ; 


= -+2s, 
2 


(13) 


Eqs. (fT^ have the same form for projection operators ttJ 
and TT^ on singlet and triplet states of isospin. Functions 
farir) in Eqs. (fT0]) - (fT^ depend on the parameters 
W^, B^, M^, where C stands for superscripts ’C’, 

’SO’ and ’T’, which are given in Refs. [Hj, [sj. In this 
work, we adjust the coupling constants in Eqs. (HI 
ini), where a and a' are indices of either singlet (s) or 
triplet (t) states. 


B. The coupled-channel equations 


and 

Nc',c{r', r) = (r', c'|r, c> (18) 

are the Hamiltonian matrix elements and the norm ma¬ 
trix elements in the channel representation, respectively. 

C. Channel states expansion in the Berggren basis 

In the present studies, any target state |ctarg) is an 
antisymmetrized state of H - 1 nucleons: 

Ictarg) = E(SDl^’'^|ctarg)|SD(^-l)) = ^ ISD^^-')) 

i i 

(19) 


Nuclear reactions can be conveniently formulated in 
the CC representation of the Schrddinger equation. The 
first step to derive the GSM-CC equations is to expand 
the GSM eigenstates in the complete basis of channel 
states {|c)} E {|cproj; Ctarg)} which contain information 
about the structure of the target and the projectile. In¬ 
dices Cproj and Ctarg denote the sets of quantum num¬ 
bers associated with the projectile and the target, re¬ 
spectively. The nuclear reaction is then described by the 
relative motion of target and projectile nuclei and the 
channel parameters, like angular momenta of the target 
and the projectile and quantum numbers of the target in¬ 
ternal excitations. In the following discussion, the heavy 
reaction participant is called a ’target’ and the light one 
a ’projectile’. Obviously, the formulation of reaction the¬ 
ory in the GSM-CC approach does not depend on this 
arbitrary choice of labels. 

The antisymmetric eigenstates of GSM-CC equations 
= = drr^ {r,c\^)\r,c) , (14) 

C 

where A is the antisymmetrization operator, 
can be expanded using the channel basis states: 


Slater determinants |SD-"^~^^) are built using a complete 
set of single-particle states of the Berggren ensemble 
which includes both resonant states and complex-energy 
scattering states. Berggren ensemble is generated by the 
single-particle potential C/basis acting on the valence nu¬ 
cleons. This ensemble is also used to generate the states 
of the projectile: 

l<(>.;cp„,> = ^(IC")®|cproj))=i(IC">®IM;j,m (20) 

where and |cproj) are radial and angular parts, re¬ 

spectively. In this expression, I is the orbital angular 
momentum of the nucleon, s its spin, j is the total angu¬ 
lar momentum and mj its projection. Hence, the basis 
state l^,Cproj) of a projectile can be written as: 

Cproj) = E l^i;cproj > (21) 

i ^ 

where Ui{r)lr = Using Eq. (HD), one can write 

the channel basis states as: 

|r,c) = E^IC^c) (22) 

• r 

I ' 

where |C^c> = i(|0r‘'>®|c)). 
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D. Hamiltonian matrix elements 

Matrix elements of the Hamiltonian Hc’^ci.r', r) and the 
norm ,r) can be derived using the expansion (1221) 

which allows to treat the antisymmetry in the projectile- 
target system. In practice, only a finite number of Slater 
determinants contribute significantly to the target state 
and thus the antisymmetry between the low-energy tar¬ 
get states and the high-energy projectile states can be 
neglected in most cases. The high-energy terms corre¬ 
spond to the channel basis states with high momentum 
k or high-i indices i > imax, where imax depends on the 
considered channel c. Hence, the expansion (1^^ can be 
splitted into low- and high-energy parts: 


Slater determinants to take into account the antisymme¬ 
try. In the second and third cases: i < imax and i' > imax 
and i > imax and i' < imax which are symmetric with re¬ 
spect to the exchange of i and i', the matrix elements are 
equal zero because Berggren states |^i;cp„j) and |</>i';cproj) 
with i > imax or i' > imax are orthogonal to all target 
states. In the last case: i > imax and i' > imax, there is 
no antisymmetry and only terms with i = i' are non-zero. 
One obtains: 


Hc'Ar'.r) 


2/i \r dr"^ 


^0 + 1) ,2 

^2 '^Ctarg 


^ „2 '^darg'Ctarg + ^C',c(r ,r) (29) 


|r,c) 



W,c) 

N 

+ E 

Uz{r) 1 


ti r 


^“^max 

r 


W-i 


N 

+ E 

u,{r) 1 




^~'^max 

r 

(23) 


where N is the number of discretized continuum states 
and imax is the index from which the antisymmetry effects 
are neglected. Equivalently, Eq. (l23l) can be written as: 


wx 1 Ut{r) 


^>= E 


\^r,c) 


2=1 


^ max 1 ?/ 

")®|C>- ^ ® Ictarg) (24) 


where fcctarg “ and the channel-channel cou¬ 

pling potential Vc\ci'f’',r) is given by: 

= UhasUr) ^ ^ (30) 

with 




_ *““”1 u^jr') Ui{r) 
E-/ ^ 


(.^2,Cproj ,Ctarg -^Ctarg ) 

(31) 


In the same way, for Nc'^df' ,r) one obtains: 


where |r) ® |c) and |<(>i;cproi) ® |ctarg) stand for non- 
antisymmetrized states. In this particular case (i > imax), 
the GSM Hamiltonian @ splits into i^proj and Htaig 
terms acting on projectile states \4‘i-,cp,oi) and target 
states |ctarg), respectively. Moreover, 

Hproj \4‘i;Cproj ) ~ ^i,Cpioj l0j;Cproi ) (25) 

.fftarg |ctarg) = ^ictarg |ctarg) • (26) 

Matrix elements of the Hamiltonian: 


Nc',c{r',r) = 


S(r - r') 


• i^targ 


+ N,,4r',r) (32) 


with: 


AT (J ) Ui{r) 

^c',C\'^ 5^) “ y. . ^c'^c\^ 1^) 

‘~T‘ . T*' nr' 


iu,(r')u,(r)^ 

J a, "darg.Ctarg ' 


ifecc(r',r) = (01' ,c'|7^|C^c> 


N 


vE 


Ui'{r') Ui{r) 
r' r 


Hc',c{i',i) 


(27) 


and of the norm: 


iVeCc(r',r) = 


N 


vE 


Ui'(r') Ui(r) 

r' r 


Nc’,c{i , i) 


(28) 


are evaluated using the expansion (l24l) . 

In the calculation of sums in Eqs. (|27|) and (l28|) . four 
cases have to be considered. In the first case: i < imax and 
i' < imax, the matrix elements are calculated in terms of 


E. Orthogonalization of the channel states 

The CC formalism leads to a generalized eigenvalue 
problem because different channel basis states are non- 
orthogonal. The non-orthogonality of channel states 
comes from the antisymmetry between the projectile and 
target states. To formulate GSM-CC equations as the 
generalized eigenvalue problem, one should express Eq. 
dT51) in the orthogonal channel basis {l^, c)^}: 

o {r, c'|r, c)^ = dc'c ■ (34) 

The transformation from the non-orthogonal channel ba¬ 
sis {|c,c)} to the orthogonal one {|r, c)^} is given by the 
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overlap operator O such that: |r,c) = 02 |r,c)^. The CC 
equations written in the orthogonal basis are: 


f. drr'^io{r',c\Ho\r,c)^-Eo{r\c\6\r,c)J 

C 

X o(r,c|^'o) = 0 , (35) 


where o{r',c'\Ho\r,c)^ = {r',c'\H\r,c), o{r',c'\d\r,c)^ = 
(r',c'|r,c), and o{'r^c\^o) = (?',cldt). The transforma¬ 
tion of this generalized eigenvalue problem into a stan¬ 
dard eigenvalue problem is achieved with a substitution: 
1$) = OlT). One obtains: 


ii: 


^2 11 J\ 


Weir) ^We'ir') 

CO 

r r' 

(r',c'|Hm|r, c). 

where 


with o{r' ,c'\H\r,c)^ = 

^ ^ ^ 

Hm = 0~^H0~^ is the modified Hamiltonian. 


Matrix elements of Hm are calculated using the expan¬ 
sion (l24l) as described in Sec. Ill PI In order to have a more 
precise treatment of the antisymmetry in the calculation 
of matrix elements of Hm, we introduce a new operator 
A: 0~2 = As-1, which is associated with the part of 0^2 
acting on the low-energy channel states. Then, instead 
of calculating the matrix elements of Hm directly, it is 
possible to calculate them as: 


f. drr'^ioir',c\H\r,c)^-Eo{r',c\r,c)^)o{r,c\^) = 0 

(36) 

with o(r,c|$) = (r,c|Oi|d>) e Wc{r)r. In the non- 
orthogonal channel basis, these CC equations become: 


Hm = (A-ri)H(A-Hi) = H-hHA-hAH-hAHA (38) 

In this formulation, the non-antisymmetrized terms are 
taken into account exactly with the identity operator. 
Inserting (|55)) in CC equations (EZD and replacing matrix 
elements (r',c'|H|r,c) using (P^)) and (I5(I1) . one obtains 
the CC equations for the reduced radial wave functions 

Wc{r)lr: 


1 ■2f‘\r 


d^(r-) l(ltl) 




^(non-loc)/ 

, ,2^c,c' VH)Wc'{r) 

rr - - - - — 




S(r - r') Wc{r) 


• i^targ 


(39) 


with the local potential vi*°'^^(r) = U\,a.sm{r) which may 
depend on the channel c, and the non-local potential: 


-e 


(non-loc)/ ! 


(r',r) = Vc',c{r',r) + (r', c'|HA|r, c> 


+ (r',c'|AH|r, c) + (r', c'|AHA|r, c) 


(40) 


The radial channel wave functions Ucir)jr are then ob¬ 
tained from the solutions of Eq. (15^ using the equation: 


source term Sdr) which would replace local t4c^°'^^(r) 
and non-local ^°'^^(r, r') potentials in (l3^ . Such an 

equivalent potential is defined by: 


V^7\r) = 

, f°° J ; T r(non-loc) / i\ / i\ 

■ L / dr \r,r )Wc'{r ) 

n' -^0 


Wc'{r) 


(42) 


and a corresponding source term is: 


Uc{r) 

r 


Wcjr) ^ y-, r 

r if Jo 


oo 


dr' r'^ (r, c\0 = AO 2 |r', c') 


Wc'jr') 

r' 

(41) 


Scij) = Ec.'{r)Y, f dr'V7?’^~^°''\r,r')wc'ir') . (43) 

c' 

Ec{r) in Eqs. (1421) . (l4^ is the smoothing function: 


F. Solution of the GSM-CC equations 

CC equations (l3^ contain a non-local potential which 
has to be treated using a generalization of the method 
of the equivalent potential [H, H^. The basic idea is 
to find the equivalent local potential the 


Ec{r) = exp -a 


Wcjr) 

Kir) 


1 - exp-a 


^asymp(^) 


Weir) 


- 1 


(44) 


to cancel divergences of the equivalent potential i^) 
close to the zeroes of wdr). In this expression: w'dr) = 
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Wc{r)lr, and is the asymptotic form of Wc{r) 

when r ~ 0. Typically, the value of a varies in the in¬ 
terval 10 < a < 100. 

With these substitutions, the GSM-CC equations (l39l) 
become: 




dr^ 


(45) 


where: = 2^{E - . Eqs. dJS]) are solved iter¬ 

atively to determine the equivalent potential, the source 
term, and the mutally orthogonal radial wave functions 
Wc(r). Starting point for solving these equations is pro¬ 
vided by a set of radial channel wave functions {uic(?')} 
obtained by the diagonalization of GSM-CC equations 
(1571) in the Bergen basis of channels. Diagonalization of 
CC equations in the Bergen basis was also considered in 


Ref. [ 5 ^. Note that it is numerically more convenient 
to express the potential r') of Eg. 1 (1391) 1 in a 

basis of harmonic oscillator states, as r') is 

short-range. For this, it is sufficient to replace all occur¬ 
rences of Berggren basis functions Ui(r)fr by harmonic 
oscillator states overlaps in Eas. (l31l33p . where 

1^1^is a harmonic oscillator state. 


III. THE RADIATIVE CAPTURE PROCESS 

We now discuss the calculation of proton/neutron ra¬ 
diative capture cross sections using the antisymmetrized 
initial and final GSM wave functions. The differential 
cross section for a proton or neutron radiative capture 
can be calculated from the Fermi golden rule, which re¬ 
lates the cross section to the matrix elements of a tran¬ 
sition operator between an initial state |i) of energy Ei 
and a final state |/) of energy Ef. The differential cross 
section is given by: 


da 1 (^'y\ ( ^^\ (\ 1 1 

dQy StT \ fc / y ?iC / \ he / 2s -I- 1 2 Jtarg + 1 


X Z 

PyTUs 


^*^V2^(2L + 1) 

L 


(?) 


L+l 


P 


L {2L + iy.\ 




J_ (^\ ( Z'l ( kuC^ \ 

27r \ fc / \ ?ic / \ he ) 


1 1 
2s -I- 1 2 Jtarg + 1 


E 

MtargiMi,, 

Pyrug 


X! ^7j ^7; ^7) (^/(^/: -^/) |•^L,ML ) 

L 


2 


(46) 


where: 

L + (47) 

In the above expressions, k^ (in units of fm”^) 

is the linear momentum of the emitted photon: 

k^ = {Ef - Ei) I {he), e^l{hc) is the electromagnetic cou¬ 

pling constant, k (in units of fm~ ) is the linear momen¬ 
tum of the incoming proton in the GM reference frame, 
(in MeV) is the reduced mass of the total system 
of A nucleons, s is the spin of the proton, Jtarg is the 
total angular momentum of the target, P = ±1 is the 
polarization of the photon, L and Mp are the multi¬ 
poles and multipole projections of the photon. More¬ 
over, D^^p{(fij,9^,0) is the Wigner D-matrix depend¬ 
ing on the angular variables 9^ and ipj of the photon. 


and is the electromagnetic transition operator. 

The final state |/) corresponds to the GSM-GG state 
|4'/(J/,M/)) of a total angular momentum J/ and a pro¬ 
jection Mf. The initial state |i) has a fixed value of the 
total angular momentum projection Mi and is denoted 
|$r(M0): 

Ji y Cg 

X (lc,,0,S,ms|(/ce,s)ice,Ws> 

^ (jce ! nT-sy -Ttarg; il^targKjce I JtB,rg)Jiy Mi) 

(48) 

where |'hi(J^,M^,Ce)) is the initial GSM-GG state with 
a total angular momentum J^ and an entrance channel 
quantum numbers Cg. Each set of quantum numbers Ce 
corresponds to a different channel c. This state can be 


= *^\/27r(2L + 1) 
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expressed in the channel basis as: 

C 


J 


Thus, the differential cross section (in units of fm^) writes: 

da 1 / ^7 'I / W \ 1 1 

dfl-y 27r \ fc / y ?ic/ \ he j 2s + 1 2 Jtarg + 1 


X Z 

P,ms, 

Mtaig,ML 


9M}^ ,p(^5 ^7? V^7; ^ 7) K'^i-^Oce ) (^Ce0s?7T-s|jcg?7ls ) (jce^s'Aarg-^targj'A-^i ) 

L Ji,Ce 


2 


X E E 

Mi,Ms, L,L', 

P,mB, Ji,J', 

Mtarg.Mi 

X ( ( J/M/|AdL',Mi,|(4-^i)c'J X (^CeOsmsIjceTOs) (4'0sms|jc;Ws> 

^ {jce^s'^targ-^targl'TiTfz) (JcgTUsdtarg-^targl^j-^i) ^ ■ (^9) 


^ 9Ml,P^^^ ^7’ ^7’ ^1^9Ml,P^^'^ ^7’ ^7’ ^ 7 ) 


J_ (El'l / E'j ( 9uC^ \ 1 _ 1 

27r \ fc / \ ftej \ ftc / 2s + 1 2 Jtarg + 1 


The operator Ml,Ml separates into an electric part 
■^L,Ml and a magnetic part m^- Formulae for the 
operators Ml are given in Appendix 171 

A. Calculation of many-body matrix elements of 
the electromagnetic operators 

The main difficulty in the calculation of matrix ele¬ 
ments comes from the infinite-range of the electromag¬ 
netic operators and the antisymmetry of the GSM-CC 
states. Indeed, direct calculation of these matrix ele¬ 
ments in the Berggren basis is not possible because they 
diverge even using the exterior complex scaling method. 

If one neglects antisymmetry in the channel state |r, c): 

|r,C> = |r>®|c> = |r>®[|Jtarg.c,Mtarg.c> ® Rc, Sc; jc, )]m 

(50) 

then the overlap between a bound state or a narrow reso¬ 
nance and a scattering state converges using the exterior 
complex-scaling method. In the above expression, Jtarg,c 
is the angular momentum of the target in a channel c 
with a projection Mtarg.c, 4 is the orbital momentum of 
the projectile, Sc its spin and jc its total angular momen¬ 
tum with a projection rrij^. The antisymmetry between 
the target and the projectile can be neglected only at 
large distances because the probability that the one-body 
state of the projectile is occupied by the target nucleon is 
the lower the smaller is the target density. In this case, 
the action of a given operator can be defined by 


considering target nucleons as distinguishable from the 
projectile nucleons: 

^Ml ~ E ^Ml (’’p (’’proj) f^proj) (51) 

i€A 

The first sum acts only on target nucleons whereas the 
second term acts on a projectile. Obviously, this approx¬ 
imation is not valid for a target in the continuum state. 

The calculation of matrix elements of the electromag¬ 
netic operators goes as follows. The matrix elements are 
expressed as the sum of a non-antisymmetrized (nas) part 
and its complement: 

(vi//||d^l|vi/,) = (vi/^lld^l|vi/,)_ 

+ ((d//||d^l|d/,)-(vi/^lld^l|vi/,)_)(52) 

The calculation of this complement is achieved by sepa¬ 
rating the operator into a short-range part and 
a long-range part . Then the symmetrized and anti¬ 
symmetrized matrix elements write: 

(vh/IIO^IIvh,) = (vh/IIO.^IIvh,) + (vh/IIO^IIvh,) (53) 

(54) 

At large distances, the antisymmetry is not crucial and 
thus the matrix element ('I'/||0> H'l'i) can be approxi¬ 
mated by (4^/||0> ||'I'i)j^^g. The remaining term is basi¬ 
cally a short-range part which can be expanded in the 










HO basis. One obtains: 

(55) 

I 


The matrix element is not antisym¬ 

metrized. We may write the operator (Eq. (IST]) 1 as: 
O^rg + *^proj! where acts only on the target state 

and Oproj O’^ the projectile state. In this case, matrix el¬ 
ements of the electromagnetic operator acting on target 
states are: 


rdrr‘^£ 


’ , , ,2Uci(r') . I ,, 

dr r ———- {r\r ) 


^ i^Cf, Scf'jjcf ) I^Ci) Sci ; jci ) ) {^Ta, ll^targll'^T'ci ) 


= J/ + 1)(2J, + 1) { 1]{JtA\6^\\Jt.) 

r' oo 

X drucf{r)uci{r) 


(56) 


where (r|r') = Sryjr^, and Ci and c/ denote initial and 
final channels, respectively. No exterior complex scaling 
is necessary to calculate the radial overlap in the above 
expression because (r) is the scattering wave function 


of a real energy and Ucf (r) is the bound state wave func¬ 
tion. Similarly, matrix elements of the electromagnetic 
operator acting on the projectile states are: 




' J , , 2 Uc,(r') 
dr r -\r|r ) 


^ {Jtcj: , I Jtc, , ) {{I 


Cf 1 ^Cf )jcA\o. 


projIK^Cij Sci)jci) 


= J/ + 1)(2J, + 1) { J 

X {Ucf,{lcf, s)jcf\\d^\\Uc,,{^Ci,s)jcA 


(57) 


r 


HO 


The antisymmetrized matrix elements ('I'/||0'^||'I'i) 
in Eq. HMD are obtained by expressing Berggren basis 
states in the HO basis. In this case, the reduced radial 
wave functions Uc(r) can be written as: 


Ucir) 


= {r\uc) 

^ Y,iAUn) {Un\Uc) 


Y,Un{r) {Un\Uc) = (J'IMc ) = 


HO\ 


(58) 


|r, c) can be expressed as: 

|r, c) = yl(|r) (8) |c)) 

= ® |c)j 

IjHO/ 'I 

= Y,{Un\r)\Un,c) = — - -\Un,c) (59) 


with |m„,c> = [ajj |Jtarg.c,Mtarg,c)l. • Hence, the 

iJCl Jc J J[^ 

CC representation of initial and final states in HO basis 
is: 


n 


(60) 


where |u„) is the radial HO state and the channel state and the antisymmetrized matrix elements of the electro- 
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magnetic operator are: 

Ci,Cf rii ,n f 

X (61) 

HO expansion is hereby justified by the fact that the 
target states are localized. 

The last many-body matrix element in Eg. (1551) : 

is calculated using Eqs. (I5S1) and (1571) 
and replacing Uc{r) by (see Eq. (l58ll l. 

IV. RESULTS OF GSM-CC CALCULATIONS 
FOR THE ^BE(P,7)®B REACTION 

GSM-CC calculations are done in COSM coordinates 
but the radiative capture cross section is expressed 
in the CM reference frame. The initial energy is: 
^(COSM) ^ ^(COSM)_^^^COSM) ^iCOSM) ^(COSM) 

and are the total energy, the projectile energy, 

and the GSM target binding energy, respectively. All 
energies are calculated in the COSM coordinate system. 
The link between the projectile energies in COSM and 
CM reference frames is given by: 

k2(ACM)\'^ 

p(COSM) ^ ^(CM) A ^ (Sroj ) A 
proj proj 2mp A-1 ^ ^ 

where is the linear momentum of the projectile. 

Energy conservation implies that the final energy is: 
^(COSM) ^ ^(COSM) ^(COSM) 

pound system binding energy in the COSM frame of ref¬ 
erence, and By = k^hc is the photon energy which does 
not depend on the chosen reference frame. 

Resonances in the spectrum of a composite A-nucleon 
system correspond to the peaks in the radiative capture 
cross section at the CM energy: Eqm = [GSM-CC]- 

E^^“^^[GSM]. Here [GSM-CC] is the GSM-CC 
energy of the resonance H’ in the nucleus A, and 
[GSM] is the GSM ground state energy of the target 
nucleus (A- 1 ). 

The cross section for a final state of the total angular 
momentum Jf is: 

r" ■ n daj^{EcM,e^,ip^) 

crjAEcM)= / d(pj / smt'^dt'^-—- 

Jo Jo dil'y 

(63) 

and the total cross section is thus: 

(y{EcM) = Jj{Egm) (64) 

•A 

In practice, one often shows the astrophysical factor: 

(65) 


which removes the exponential dependence of the cross 
section at low energies due to the Coulomb barrier, r] in 
(1651) is the Sommerfeld parameter: r] = {mZiZ 2 )/{h'^k), 
where Zi and Z 2 are the proton numbers of the projectile 
and target nuclei. 

A. Parameters of CSM calculations in ^Be and ®B 

The model space in ^Be and ®B is limited by the core of 
^He. The core is described by a WS potential (see Table 
HI for each considered partial wave: / = 0,1 and 2. The 
radius of the Coulomb potential is Tc = 2.8 fm. To deter- 


Parameter 

Protons 

Neutrons 

a 

0.65 fm 

0.65 fm 

Ro 

2.0 fm 

2.0 fm 

Vo{l = 0) 

61.5 MeV 

70.6735 MeV 

Uso(l = 0) 

0 MeV 

0 MeV 

VoQ = 1) 

44.3967 MeV 

70.6734 MeV 

Ubo(1 = 1) 

7.80188 MeV 

7.86276 MeV 

VAl = 2) 

44.3967 MeV 

0 MeV 

Uso(l = 2) 

7.80188 MeV 

0 MeV 


TABLE 1. Parameters of the WS potential of ^He core used 
in the GSM and GSM-CC description of ^Be and ®B. 

mine Berggren ensemble, one calculates first the single¬ 
particle bound and resonance states of the basis generat¬ 
ing WS potential for all chosen partial waves Then, 

for each one selects the contour in a fourth 

quadrant of the complex fc-plane. All (/,_))-scattering 
states in this ensemble belong to . The precise form 
of the contour is unimportant providing that all selected 
single-particle resonances for a given (Z,j) lie between 
this contour and the real fc-axis for 7?.(fc) > 0. For each 
(/,j), the set of all resonant states and scattering states 
on forms a complete single-particle basis. 

In the present case, valence nucleons can occupy the 
0 p 3/2 and Opi /2 discrete single-particle states and sev¬ 
eral non-resonant single-particle continuum states on 
discretized contours: Et , , Et , Et , E\ and 
Each contour consists of three segments join¬ 
ing the points: fcmin = 0.0, fcpeak = 0.15 - f0.14 fm”^, 
^middle = 0.3 fm“^ and /cmax = 2.0 fm~^, and each segment 
is discretized with 10 points. Hence, GSM and GSM-CC 
calculations are done in 152 shells: 31 P 3/2 and pi /2 shells, 
and 30 Si/ 2 , ^ 3/2 and ^ 5/2 shells. The GSM basis is trun¬ 
cated so as to reduce the size of the GSM Hamiltonian 
matrix. For this, the occupation of P 3/2 and P 2/2 scat¬ 
tering states in basis Slater determinants is limited to 
two particles, while the occupation of S 1 / 2 , ^ 5/2 and ^ 3/2 
scattering states is limited to one particle only. The lat¬ 
ter truncation is justified by the fact that GSM target 
states virtually only consists of P 3/2 and pi /2 states, S 1 / 2 , 
^ 5/2 and ^ 3/2 states occurring only in the partial wave 
decomposition of the proton or neutron projectile. 

Parameters of the Hamiltonian, which were adjusted to 


S{Ecm) = cr{EcM)EcMe^^^ 
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reproduce binding energies of low-lying states in ^Be and 
®B, are given in TableHIl In GSM calculations, the ground 


Parameter 

Value [MeV] 

Kt 

V°t 

yC 

*^t,s 

Vt®t° 

Vt^t 

4.00906 

-3.22579 

2.22077 

-9.51008 

-1448.32 

0 

15.3946 

-15.4834 


TABLE II. Parameters of the FHT interaction in GSM and 
GSM-CC calculations in ^Be and ®B. The superscripts “C”, 
“SO” and “T” stand for ’central’, ’spin-orbit’ and ’tensor’, 
respectively, and the indices “s” and “t” stand for ’singlet’ 
and ’triplet’. 

state of ^Be is bound with respect to '^He by 9.378 MeV, 
close to the experimental value i?exp = 9.304 MeV. Re¬ 
action channels in GSM-CC calculations are obtained by 
the coupling of the ground state 3/2^ and the first ex¬ 
cited state 1/2” of ^Be with the proton partial waves: 

■51/2) Pi/2j P3/2j <^3/2 ^^d ^ 5 / 2 - 

Discrete states of a composite system ®B are 2\ bound 
state, and I^, 3^, I 2 resonances. Missing reaction chan¬ 
nels in GSM-CC lead to a small difference between GSM 
and GSM-CC energies for these states. To correct this 
deficiency, the channel-channel coupling potentials 14,c' 
in GSM-CC have been adjusted for each considered state 
of ®B. The new potentials are: 14, c' = c( J’’’)I4,c') where 
the multiplicative corrective factors are: c(2j^) = 1.0133, 
c(li;) = 1.0602, and c(3t) = 1.0233. 

B. The astrophysical 5-factor for ^Be(p, 7 )®B 
reaction 

The description of electromagnetic transitions requires 
effective charges for proton and neutron. For El transi¬ 
tions, the standard values are [^ : 

eeff=e(l-|) ; = 

where Z and A are the proton number and the total 
number of nucleons, respectively. The standard values 
for E2 transitions are: 

There are no effective charges for Ml transitions. In the 
present work, we use these standard values for El and 
E2 effective charges. One should keep in mind however, 
that the effective charges extracted experimentally show 
often significant deviations from the standard values . 

Proton separation energy in the ground state of ^Be 
is S-p = 5.6 MeV. The final nucleus ®B has one weakly 



Ecu (MeV) 

FIG. 1. Plot of the El astrophysical factor for the ^Be(p, 7 )®B 
reaction. The solid line represents the exact, fully antisym¬ 
metrized GSM-CG calculation with both the ground state 
= 3/2i and the first excited state J” = l/2i of ^Be tar¬ 
get included. The dashed line shows results of the GSM-GC 
calculations if the first excited state of ^Be is omitted. For 
more details, see the description in the text. 



^CM (MeV) 

FIG. 2. The same as in Fig[T]but for the Ml transitions. The 
two peaks correspond to the l/ and 3i resonances of ®B. 

bound state 2\ below the proton emission threshold. 
Experimental proton separation energy in this state 
Sp = 0.1375 MeV agrees well with the calculated value 
= 0.137 MeV. The 1); and 3^ resonance peaks 
should be seen in Ml transitions. The 1| resonance could 
also be seen in E2 transitions. 

All relevant El, Ml, E2 transitions from the initial 
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FIG. 3. The same as in Fig[T|but for the E2 transitions. The 
two peaks correspond to the It and Ij resonances of ®B. 



FIG. 4. Plot of the total astrophysical factor for the 
^Be(p, 7 )®B reaction. Data are taken from Refs.[^ and [^ . 
The solid line represents the exact, fully antisymmetrized 
GSM-CG calculation including both the ground state = 
3/2i and the first excited state = l/2i of ^Be target. Gal- 
culations neglecting the first excited state of the target are 
shown with the dashed line. For more details see the descrip¬ 
tion in the text. 


continuum states {Ji = 1^,2^,3‘^) in ®B to the final 
bound state Jf = 2^ state have been included. Figs. [T][3] 
show the separate contributions to the total ^-factor in 
^Be(p, 7 )®B reaction: for El transitions (Fig. [IJ, 

for Ml transitions (Fig. [2), and for E2 transi¬ 
tions. The solid lines in Figs. [T][3]show results of the fully 


antisymmetrized GSM-CC calculations with both ground 
and first excited states of the ’’Be target included. The 
dashed line in these figures correspond to GSM-CC cal¬ 
culations neglecting the 1/2^ first excited state in ’Be. 

There is no resonant contribution in El transitions. In¬ 
cluding the first excited state of the target lowers 5^’ by 
less than - 5% for Ecu < 2.5 MeV. On the contrary, the 
Ml contribution to the ^-factor increases significantly in 
the region of 1^ resonance if the the excited state of the 
target is included (see Fig. n. One can see 1^ and 3^ res¬ 
onances of ®B at Ecm = 0.79 MeV and Ecu = 2.34 MeV, 
respectively. These resonances are observed experimen¬ 
tally at Ecm = 0.632 MeV and Ecm = 2.182 MeV, respec¬ 
tively. The E2 transitions contribute little to the S- 
factor. is ~10~^ smaller than 5^’ and and in¬ 
creases by less than ~ 10% for CM energies in the region 
of !(■ and IJ resonances. The resonance I 2 has not yet 
been seen experimentally. 

The calculated total S'-factor is compared with the ex¬ 
perimental data [U [ 2 ^ in Fig. H) Below Ecm = 1 MeV, 
the agreement with the data is good if both the ground 
state of ’Be and its first excited state are included. The 
value of the E-factor at zero energy, is 

23.214 b-eV and the slope, dS/dEcM\EcM=o^ is 37.921b. 
The accepted experimental value of the E-factor is 
20.9±0.6 b-eV, slightly below the GSM-CC results. 

At higher energies, GSM-CC results overshoot the ex¬ 
perimental data. This feature could be due to the ab¬ 
sence of higher lying discrete and continuum states of 
’Be target in the channel basis. Indeed, in the present 
case, GSM and GSM-CC calculations with uncorrected 
channel-channel coupling potentials 14,c' do not give the 
same spectra and binding energies of ’Be and ®B and the 
small multiplicative correction factors are necessary. 

The long wavelength approximation simplifies the cal¬ 
culation of matrix elements of the electromagnetic tran¬ 
sitions. The quality of this approximation and the 
role of the antisymmetry of initial and final states is 
tested in Figs. [5] and [6l Only the ground state of 
’Be is taken into account. To correct GSM-CC calcu¬ 
lations for the missing channels in this case, the channel- 
channel coupling potentials 14,c' have been slightly cor¬ 
rected: 14 ,c' = c( J’^)I 4 ,c '5 and the multiplicative cor¬ 
rective factors are: c{2\) = 1.0122, c(l|) = 1.0668, and 
c(3t) = 1.0225. 

At low energies (Ecm < 1-5 MeV), both the long wave¬ 
length approximation and the antisymmetrization in the 
calculation of El transition matrix elements does not 
change results significantly (see Fig. [5]). Both approx¬ 
imations become worse at higher energies but even at 
Ecm = 2.5 MeV the error is only ~ 10%. 

The astrophysical factor for Ml transitions is shown 
in Fig. [6l The antisymmetrization of the initial and 
final states lowers the value of E*^’ by a factor ~ 2 at the 
resonance peaks. The long wavelength approximation 
does not change E'^’. 
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Ecu (MeV) 

FIG. 5. Plot of the El astrophysical factor for the ^Be(p, 7 )®B 
reaction. The solid line represents the exact, fully an¬ 
tisymmetrized GSM-CC calculation. The calculations in 
the long wavelength approximation are represented by the 
dashed and dotted lines in the fully antisymmetrized and non- 
antisymmetrized cases, respectively. For more details, see the 
description in the text. 



^CM (MeV) 

FIG. 6. The same as in FigObut for the Ml transitions. The 
two peaks correspond to the It and 3t resonances of ®B. 


V. RESULTS OF GSM-CC CALCULATIONS 
FOR THE ^LI(N,7)®LI REACTION 

A. Parameters of GSM calculations in ^Li and ®Li 

^Li(n, 7 )®Li is the mirror reaction of ^Be(p, 7 )®B and 
will be described in the same model space.The WS po¬ 


tential of "^He core is given in Table IIIII The radius 
of the Coulomb potential is = 2.8 fm. Valence nu- 


Parameter 

Protons 

Neutrons 

a 

0.65 fm 

0.65 fm 

Ro 

2.0 fm 

2.0 fm 

V ^{1 = 0) 

71.0752 MeV 

43.6438 MeV 

Vso(? = 0) 

0 MeV 

0 MeV 

Vo{l = 1) 

71.0752 MeV 

43.6438 MeV 

= 1) 

7.90622 MeV 

7.84517 MeV 

V ^{1 = 2) 

0 MeV 

43.6438 MeV 

Vsoil = 2) 

0 MeV 

0 MeV 


TABLE III. Parameters of the WS potential of “^He core used 
in the GSM and GSM-CG description of ^Li and ®Li nuclei. 

cleons occupy OP 3/2 and Opi /2 discrete single-particle 
states and non-resonant single-particle continuum states 
on discretized contours: 1 and 

Each contour consists of three segments join¬ 
ing the points: fcmin = 0.0, fcpeak = 0.15 -10.14 fm”^, 
^middle = 0.3 fm”^ and fcmax = 2.0 fm~^, and each segment 
is discretized by 10 points. 

Parameters of the FHT Hamiltonian in ^Li and ®Li 
are given in Table HVl In GSM, the ground state of ^Li 


Parameter 

Value [MeV] 


4.03185 


-4.95286 


2.23361 


-7.63465 


-1456.55 

vi!? 

0 

Vt^t 

15.4822 

VI3 

-15.5716 


TABLE IV. Parameters of the FHT interaction for GSM and 
GSM-CG calculations in ^Li and ®Li. For more details, see 
the caption of Table im 

is bound by 11.228 MeV with respect to ^He, i.e. close 
to the experimental value (Eexp = 10.948 MeV). Reac¬ 
tion channels are obtained by the coupling of the ground 
state 3/2 and the first excited state 1/2 of ^Li with the 
proton partial waves: Si/ 2 , Pi/ 2 , P 3 / 2 , ^ 3/2 and ^ 5 / 2 . 

Discrete states of a composite system ®Li are 2)', 
Ij bound states, and 3^ resonance. To correct for 
missing reaction channels in GSM-CC calculations, the 
channel-channel coupling potentials I4,c' have been mod¬ 
ified and new potentials are: V/.c' = c(J’^)Vc,c 7 with 
c(2+) = 1.03705, c(lt) = 1.04805 and c(3t) = 1.03205. 

B. ^Li(n, 7 )®Li cross section 

Neutron separation energy in the ground state of 
^Li is Sn = 7.25 MeV. The final nucleus ®Li has two 
bound states J'^ = 2 ^ and 1^ below the neutron emission 
threshold. Neutron separation energy from the ground 
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state and excited states are Sn = 2.03262 MeV and 
S^= ... MeV, respectively. The calculated neutron sepa¬ 
ration energies in these two states are = 2.032 MeV 
and = ... MeV, in good agreement with the experi¬ 
mental data. The 3^ resonance peak can be seen in Ml 
and E2 transitions. 



Ecu (MeV) 

FIG. 7. The same as in Fig[T]but for the reaction ^Li(p, 7 )®Li 
reaction. For more details, see the description in the text. 



(MeV) 

FIG. 8. The same as in Fig[T] but for Ml transitions in 
^Li(p, 7 )®Li reaction. The peak corresponds to the 3i res¬ 
onance in ®Li. 

Figs. SI show the El, Ml and E2 cross sections 
for ^Li(n, 7 )®Li reaction. The solid lines in Figs. S El 
show results of the fully antisymmetrized GSM-CC cal¬ 
culations with both ground and first excited states of the 


Xl0“^ 



Ecu (MeV) 

FIG. 9. The same as in FigObut for the E2 transitions. The 
two peaks correspond to the 3i and ij resonances of ®Li. 



FIG. 10. Plot of the total cross section for the ^Li(n, 7 )®Li 
reaction. Data are taken from Ref. [^. The solid line repre¬ 
sents the exact, fully antisymmetrized GSM-CG calculation. 
Calculations in the long wavelength approximation are repre¬ 
sented by the dashed and dotted lines in the antisymmetrized 
and non-antisymmetrized cases, respectively. For more details 
see the description in the text. 

^Li included. The dashed line in these figures correspond 
to GSM-CC calculations neglecting the 1/2“ first excited 
state in ^Li. Including the first excited state of the target 
lowers El contribution to the neutron radiative capture 
cross-section by ~ 20% for Ecm < 1 MeV. 

The Ml contribution to the cross-section increases by 
~ 25% in the region of 3^ resonance if the excited state 
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of the target is included (see Fig. |8l). One can see that 
the calculated 3^ resonance is at the experimental value 
of energy: Ecm = 0.223 MeV. 

E2 transitions contribute very little to the neutron ra¬ 
diative capture cross section. The E2 contribution is 
three orders of magnitude smaller than El and Ml con¬ 
tributions. The role of the excited state of the target is 
very important. It increases the contribution from E2 
transitions by a factor -- 3 in the region of 3^ resonance. 
At the I 2 resonance, the excited state enhances the E2 
contribution by about one order of magnitude. The cal¬ 
culated energy of this resonance is lower than seen ex¬ 
perimentally. 

The total neutron radiative capture cross section is 
compared with the experimental data in Fig. [101 
GSM-CC calculation underestimates the data of Imhof 
et al [ 5 ^. The extrapolation of the calculated neutron 
radiative capture cross section at low Ecm is done using 
the expansion: 

a(EcM) = 4^-2.360 + 3.387/^ (68) 

n/Ecm 

which yields: (j(GSM-cc) _ 25.41 /^barn at Ecm = 25 keV. 



Ecu (MeV) 

FIG. 11. The same as in Fig|5]but for ^Li(p, 7 )®Li reaction. 
For more details, see the description in the text. 

The long wavelength approximation and the role of 
the antisymmetry of initial and final states in the calcu¬ 
lation of matrix elements of the electromagnetic transi¬ 
tions is tested in Figs. [TTl [TH Only the ground state 
of ^Li is taken into account in these tests. To correct 
GSM-CC calculations for the missing channels in this 
case, the channel-channel coupling potentials 14 ,c' have 
been corrected: 14,c' = c( J’^)14,c', and the multiplicative 
corrective factors are: c{2\) = 1.038, c(l(') = 1.0594, and 
c(3(;) = 1.032. 

At low energies (Ecm < 1-2 MeV), both the long wave¬ 
length approximation and the antisymmetry of initial 



FIG. 12. The same as in Fig llll but for the Ml transitions. 
The peak corresponds to the 31 resonance in ®Li. 


and final states in the calculation of El transition ma¬ 
trix elements does not change results significantly (see 
Fig. fTlT) . Also Ml transition matrix elements are in¬ 
sensible to the long wavelength approximation (see Fig. 
dsi). On the contrary, the antisymmetrization is essential, 
decreasing the Ml contribution to the neutron radiative 
capture cross section by a factor ~ 4 in the region of 3(^ 
resonance. 


VI. CONCLUSIONS 

The GSM in the coupled-channel representation opens 
a possibility for the unified description of low-energy 
nuclear structure and reactions using the same Hamil¬ 
tonian. While both GSM and GSM-CC can describe 
energies, widths and wave functions of the many-body 
states, the GSM-CC can in addition yield reaction cross- 
sections. Combined application of GSM and GSM-CC 
to describe energies of resonant states allows to test the 
exactitude of calculated cross-sections for a given many- 
body Hamiltonian. 

In this work, we have presented in details the GSM 
in the coupled channel representation and applied it for 
the description of the low-energy proton and neutron ra¬ 
diative capture processes on mirror targets ^Be and ^Li, 
respectively. The interaction between valence nucleons 
in this calculation was modelled by the finite-range two- 
body FHT interaction. 

The convergence of GSM-CC calculations has been 
checked by comparing GSM and GSM-CC results for ®B 
and ®Li states. In a given single-particle model space, the 
GSM-CC calculation with the reaction channels which 
are constructed using selected many-body states of the 
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target nucleus (^Be or ^Li in our case), can be considered 
reliable if the GSM-CC eigenvalues for a combined sys¬ 
tem (®B or ®Li in our case) approximate well results of a 
direct diagonalization of the GSM Hamiltonian matrix in 
the same single-particle model space. In such a case, the 
configuration mixing in GSM-CC and GSM wave func¬ 
tions are equivalent and one does not need to include 
additional states of the target nucleus to reach the many- 
body completeness in GSM-CC calculation. Only in this 
case, the unified description of nuclear structure and re¬ 
actions with the same many-body Hamiltonian and the 
same model space is reached. In the studied case, the 
GSM and GSM-CC spectra were close but not identi¬ 
cal so the small renormalization of the channel-channel 
coupling potentials was necessary to compensate for the 
missing channels made of the higher lying discrete and/or 
continuum states of the target. 

There are two important aspects in this GSM-CC cal¬ 
culations which have been studied carefully. The first 
one is the antisymmetry of initial and final states in the 
calculation of matrix elements of the electromagnetic op¬ 
erators. It was found that the antisymmetry is crucial 
in Ml transitions in the region of resonances. At ener¬ 
gies of astrophysical importance, the error introduced by 
neglecting the antisymmetry is however small. The sec¬ 
ond aspect is the role of the excited state 1/2“ of the 
target. At Acm ~ 0, the radiative capture cross sections 


in ’^Be(p, 7 )®B and ^Li(n, 7 )®Li reactions are slightly im¬ 
pacted by the excited state of target. However, in the 
region of resonances and at higher energies the excited 
l/2j state in ^Be and ’’Li turns out to be crucial. As 
compared to '’Be(p, 7 )®B reaction, the '’Li(n, 7 )®Li reac¬ 
tion is less sensitive to the first excited state of the tar¬ 
get but more sensitive to the antisymmetry of initial and 
final states in the calculation of matrix elements of elec¬ 
tromagnetic transition operators. The long wavelength 
approximation in the transition matrix elements changes 
mainly the El contribution to the radiative capture cross 
section. 


Appendix A: Matrix elements of the 
electromagnetic transition operators 

The matrix elements related to electric and magnetic 
transitions will be considered with and without long- 
wavelength approximation. The operators involving the 
exact and approximate electromagnetic field can be found 
in Ref. (60|, whose matrix elements can be derived 
straightforwardly from the Wigner-Eckhart theorem and 
standard manipulations of gradients of spherical har¬ 
monics coupled to angular momenta [6l|. The opera¬ 
tor Ml,Ml separates into electric and magnetic 
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where i runs over all considered nucleons and li and Si are the above expression, Y^^{Vt) is a spherical harmonics, 
the orbital and spin angular momenta, respectively. In 
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S'l is the Ricatti-Bessel function, ri, fli are radial and an¬ 
gular coordinates of the nucleon i, and Un = filvi. More¬ 
over, Ci is the dimensionless charge of the nucleon i (e^ = 1 
for a proton and 0 for a neutron), gl is the dimensionless 
magnetic spin moment of the nucleon i = 5.5857 for 
a proton and -3.8263 for a neutron), mp<? (in units of 
MeV) is the mass of the proton, and g\ is the dimension¬ 
less magnetic orbital momentum of the nucleon i times 


L + 1 {g\ = 2 for a proton and 0 for a neutron). 

In the long wavelength approximation, the expressions 
dAH) and ()A2p become: 

(A3) 


AbM _ sr 


L + 1 




he 


2mpC^ 


he 


2mpC^ 


^ VA(2L + l)rf-i ^ [Y^-^(fl.) ® U 
^ VL(2L + 


(A4) 


Eqs. (lUll-dAll have been written in such a way that only 
one-body operators appear in each summation. Matrix 
elements of these operators are calculated in a standard 
way using the Wigner-Eckhart theorem. 
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